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Massive black hole binaries are key targets for the space based gravitational wave interferometer 
LISA. Several studies have investigated how LISA observations could be used to constrain the 
parameters of these systems. Until recently, most of these studies have ignored the higher harmonic 
corrections to the waveforms. Here we analyze the effects of the higher harmonics in more detail 
by performing extensive Monte Carlo simulations. We pay particular attention to how the higher 
harmonics impact parameter correlations, and show that the additional harmonics help mitigate 
the impact of having two laser links fail, by allowing for an instantaneous measurement of the 
gravitational wave polarization with a single interferometer channel. By looking at parameter 
correlations we are able to explain why certain mass ratios provide dramatic improvements in 
certain parameter estimations, and illustrate how the improved polarization measurement improves 
the prospects for single interferometer operation. 

PACS numbers:04.25.Nx, 04.30.Db, 04.80. Cc 

I. INTRODUCTION 

Massive black hole binaries are expected to be some of the brightest gravitational wave (GW) sources for the 
planned Laser Interferometer Space Antenna (LISA) This joint ESA-NASA mission will search for GWs is the 
band 10~ 5 < //Hz < 1. While the event rates are uncertain, it is likely that there will be at least a few events per 
year. Massive black hole binaries are a very promising source of GWs as we expect them to have integrated signal 
to noise ratios (SNRs) in the hundreds to thousands. As well as being very bright, there should be little confusion 
with other sources (a problem which plagues galactic binary extraction, and to a lesser extent the extreme mass ratio 
inspiral sources). 

There have been many studies looking at parameter estimation for spinning and non-spinning binaries 0, H, 0, IH 
El 0j B i i_ 13] j and, more recently, on the development of search algorithms [ill, EH EU 0, EH • With the exception 
of Refs. 

these studies have focused on the contribution from the dominant second harmonic of the orbital 
frequency, and have neglected the influence of the additional harmonics that appear at higher post-Newtonian order. 
The studies that did consider the impact of these higher harmonic corrections (HHCs) saw improvements in parameter 
estimation, but o nly a few cases were considered. More extensive studies of the impact of HHCs have appeared in the 
last year [l(| E3,EIj], and it has become apparent the effects can vary significantly from source to source. 

Our goal here is to consider a wider sample of systems, and to elucidate the mechanism by which the HHCs affect 
parameter estimation. Our most significant finding is that HHCs greatly mitigate the effects of any hardware failures 
that reduce the instrument to a single interferometry channel [19(. Indeed, when HHCs are included, the science 
performance with a single channel is comparable to what was found for two channels when HHCs are neglected. 

We performed extensive Monte Carlo simulation for various redshifted chirp masses, with coalescence time t c = 1, 
lumiosity distance Dr, = lOGpc and all other parameters chosen at random. Results for other distances can be 
obtained by multiplying our parameter uncertainities by lOGpc/ Dr,. As sky resolution is an important quantity, we 
considered 20,000 different sky locations per chirp mass and per mass ratio. Finally, for each chirp mass, we run 
separate simulations for mass ratios of 1 and 10. 



A. Outline of the paper 

The paper is structured as follows. In Sec ([H]) we outline the form of the gravitational waveform at the detector 
with higher harmonics corrections. We also give a brief description of the low frequency approximation Q for the 
LISA response. Sec. contains a brief outline of the main tools used for parameter estimation. We finish this 
section with an outline of the Monte Carlo simulations we carried out. In Sec pV[) we present results of our analysis 
on the effect of the higher harmonic corrections for parameter estimation using different redshifted chirp masses and 
mass ratios. This is followed in Sec (|V)) by an investigation into what correlation breaking between parameters is 
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responsible for the increase in parameter estimation. The final main section, Sec (|VI[) , deals with the effect of losing 
interferometry links during the mission. 

II. THE GRAVITATIONAL WAVEFORM 

In the low frequency approximation, we can think of LISA as being composed of two orthogonal 90° detectors. The 
strain of the gravitational wave (GW) in each channel of the LISA detector with both polarizations is given by 

h(t) = h + (at))p + + hAm)F x , (i) 

where the phase shifted time parameter is 

f(f) = t - i? ffi sin cos (a(i) - 0) . (2) 

Here, Rq = IAU « 500 sees is the radial distance to the detector guiding center, (0, <j>) are the position angles of the 
source in the sky, a(t) = 2-Kf m t + k, f m = 1/year is the LISA modulation frequency and k gives the initial ecliptic 
longitude of the guiding center. The GW polarizations up to 2-PN order in amplitude corrections are defined by [201 ] 
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Here m = mi + mi is the total mass of the binary, r\ — m\mijm 2 is the reduced mass ratio and Dl is the luminosity 
distance of the source. The invariant PN velocity parameter is defined by, x = (Gmw/c 3 ) 2 3 , where u> = d<$>o r b/dt is 
the 2 PN order orbital frequency for a circular orbit and $ or b = tp° rb — 4> rb(t) is the orbital phase which is defined as 

* / \ nrh 1 f^/« /3715 55 \ 3?r 1/4 / 9275495 284875 1855 2 \ ^ 1/R 1 

$ orb (t) = ip°; b - - 1 e j/s + h —i] e 3/8 e 1/4 + 1 n + — n 2 e 1/8 y , 4 

v ' Vc r)\ V8064 96 7 4 V 14 450688 258048 ' 2048 ' / J ' v ; 

where the quantity O is related to the time to coalescence of the wave, i c , by 

e(*) = ^(*o-t). (5) 

and is the orbital phase of the wave at coalescence. All GW phases are then twice the orbital value. For the rest 
of the paper, we will work with GW phases. 

In Eqn ©, the functions fll x contain the PN corrections to the amplitude and the extra phase harmonics. The 

restricted PN waveform corresponds to keeping just the H+ x terms. We should note here that contained in the half 

integer H^ x terms is a factor 5m — mi — mi. This term has the effect of killing all the odd phase harmonics in the 
equal mass case. We can also see the extra frequency harmonics arising due to the x n terms. 

Using the WMAP values of (fl R , Q M , fi A ) = (4.9 x 10" 5 , 0.27, 0.73) and a Hubble's constant of H =71 km/s/Mpc, 
the relation between redshift, z, and luminosity distance, Dl-, is given by 



. . nK(i + /) 4 + o M (i + 2') 3 + o A l 1/2 . (6) 

-no Jo L J 

The functions F + ' x are the beam pattern functions of the detector given in the low frequency approximation by 
F+(t; xP, 0, 0, A) = \ [cos(2^)^+(i; 0, 0, A) - sin(2^)£» x (f; 0, 0, A)] , (7) 

^ x (t; V, 0, 0, A) = 1 [sm(2iP)D+(t; 0, 0, A) + cos(2^)D x (t; 0, </>, A)] , (8) 

where ip is the polarization angle of the wave and A = or 7r/4 defines the two-arm combination of LISA from which 
the strain is coming. The detector pattern functions are given by [2ll | 

D+{t) = ^ - 36sin 2 (0)sin(2a(i) - 2A) + (3 + cos(20)) ( cos(2<£) { 9sin(2A) - sin(4a(t) - 2A)| (9) 
+ sin(2c/>) | cos(4a(t) - 2A) - 9cos(2A)| \ - 4V3sin(20) ( sin(3a(f) - 2A-0) 



3sin(a(i) -2\ + <j>)\ 
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D x (t) = — \ V3 cos(6») f 9 cos(2A - 20) - cos(4a(t) - 2A - 20) 
16 L V 

-6 sin(6>) ( cos(3a(t) - 2A - 0) + 3 cos(a(i) - 2A + 0)) 



(10) 



For two Schwarzschild black holes, the above equations governing the evolution of the phase break down even before 
we reach the last stable circular orbit (LSO) at R = 6M. Because of this, we terminate the waveforms at R = 7M. 

The low frequency approximation is an extremely good fit to the full detector response at frequencies of < 3 
mHz [21]. In this case the two channel formalism originally derived by Cutler corresponds to the construction of 
optimal orthogonal time delay interferometry (TDI) variables {^4, E} using the unequal-arm Michelson TDI variables 
{X, Y, Z} according to 



A = X , E= (X + 2Y)/V3. 
In later sections we will refer to the one channel case as X and the two channel as AE. 



(11) 



III. ESTIMATING PARAMETER ERRORS USING THE FISHER MATRIX. 



One of the main tools used in the GW community for the estimation of parameter errors is the Fisher information 
matrix (FIM). In the high SNR limit, the inverse of the FIM gives the variance-covariance matrix. The square root 
of the diagonal elements of the inverse FIM give a 1-tr estimation of the error in our parameter estimation. For 
LISA, it has been shown in a number of cases that the FIM is a good parameter error estimator when compared 
with Markov Chain Monte Carlo techniques. For this particular problem, the parameter set we chose to work with is 
A M = {ln(M c ), ln(/x),ln(t c ), cos 6,<f),ln(D L ), cost, tpd 3 } where M c — mrf^ is the chirp mass, [i = mr] is the reduced 
mass and all other parameters have been previously defined. It has been customary to divide this parameter set into 
intrinsic (parameters that effect the dynamics of the system) and extrinsic parameters (those that are more relative 
to the detector). In the parameter set defined above, we consider the final four parameters to be extrinsic. The Fisher 
matrix is thus defined by 



r = 



dh 
dA 1 



dh \ 
9AT/ ' 
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where h = h(f) is an un-normalizcd template and the angular brackets denote the inner product 

df 



(h\a) 



Sn(f) 



h(fys*{f) + h*(f)s(f) 



(13) 



where a tilde denotes a Fourier transform and an asterisk denotes a complex conjugate. The quantity S n (f) — 
ST str (f) + Sf° nf {f) is the one-sided noise spectral density of the detector, which is a combination of instrumental 
and galactic confusion noise. For the instrumental noise we use the expression given by [22j 
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where L = 5 x 10 6 km is the arm-length for LISA, S£ os (/) = 4 x 10" 



2 1 Hz and S° cc {f) = 9 x lO" 30 m 2 /s i /Hz 



are the position and acceleration noise respectively. The quantity /* = l/(27rL) is the mean transfer frequency for 
the LISA arm. Notice that the final term in the expression has the effect of reddening the noise below 10~ 4 Hz to 
account for the fact that we may not be able to achieve the desired noise spectral density as we approach 10~ 5 Hz. 
For the galactic confusion we use the following confusion noise estimate derived from a Nelemans, Yungelson, Zwart 
(NYZ) galactic foreground model [13, H3] 
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where the confusion noise has units of Hz^ 1 . When two channels are available the total FIM is the sum of the FIM 
in each channel, i.e. = Tjj + Tjj . 

For this study the Fisher matrix is calculated numerically, where the derivatives of the waveforms are calculated 
using the central difference equation 

dh _ h(X v + AA") - h(\" - AX") 

dX^^ 2AA" ' (1 ' 

for all parameters other than M c and [i in the equal mass case. The central difference equation in un-applicable for 
these two parameters when we have equal masses, as it is a degenerate space. Certain shifts in one parameter, while 
holding the other constant, leads to un-astrophysical individual masses. Therefore, in the equal mass case, for M c 
and \i we use 

dh _ h{M c + AM C ) - h{M c ) 

> V-L ') 



dM c AM C 
and 

dh _ h(fi) — h(fi — A/it) 

<9/U AjJL 



(18) 



For all other cases, we revert to the central difference formula. 

One of the quantities that we are most interested in investigating is the error box on the sky for a particular source. 
We therefore define the positional error on the sky as 



Afl = 2nyJ S ee E^ - (S^) 2 , (19) 
where the elements of the variance-covariance matrix are given by 

= (AA*AA J ) = (r^)" 1 , (20) 
and the three main quantities in the sky resolution expression are given by 

S ee = (A cos 6»A cos 6») , (21) 



£^ = {A(j)A(j)) , (22) 

and 

£°* = <Acos0A0). (23) 

Finally, we also quote the optimal signal to noise ratio (SNR) for each source. In each individual detector this is 
defined by 

A = (h i \h i ). (24) 
When we use both LISA detectors the total SNR is given by 

9 = \fp 2 i + P 2 m (25) 

For the Monte Carlo simulation, we used 2 x 10 4 points, where the redshifted chirp mass, mass ratio, time to 
coalescence and luminosity distance were held constant and all other parameters were varied. To investigate as many 
cases as possible, we chose chirp masses of M c (z) — {10 8 , 10 7 , 10 6 , 10 5 , 10 4 } M Q . For each chirp mass, we ran a Monte 
Carlo for mass ratios of 1 and 10. In all cases, the sources were put at a constant distance of 10 Gpc. We assume the 
time of observation is 1 year, and the time to coalescence in each case is 0.999 years. To be consistent, we chose not 
to evolve our templates beyond 5 mHz as we know that we can not trust the low frequency approximation beyond 
this value. For cases where coalescence is seen (i.e. the coalescence frequency of the maximum harmonic is less than 
5 mHz) the templates were terminated once the distance between the bodies reached 7M. 
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FIG. 1: A comparison of parameter extraction for equal mass binaries using restricted (solid line) and higher harmonic corrected 
(dashed line) waveforms, as a function of redshifted chirp mass and at a distance of 10 Gpc. The values quoted are median 
values from a 2 x 10 4 point Monte Carlo simulation. 




MJz) 

FIG. 2: A comparison of parameter extraction for binaries with a mass ratio of 10 using restricted (solid line) and higher 
harmonic corrected (dashed line) waveforms, as a function of redshifted chirp mass and at a distance of 10 Gpc. The values 
quoted are median values from a 2 x 10 4 point Monte Carlo simulation. 
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IV. THE EFFECT OF HIGHER HARMONIC CORRECTIONS ON PARAMETER ESTIMATION. 

For clarity, we will treat the equal and unequal mass results separately. In both cases, we focus on the most 
interesting parameters from an astronomical point of view. We focus on the two mass parameters, the time of 
coalescence, the luminosity distance and the sky resolution. For completeness, we also include information on the 
signal to noise ratios. We should also point out that due to the large tails in the distributions, median values are 
more informative than mean values. 

A. Equal mass binaries. 

In Fig {T]) we plot the median results from the Monte Carlo simulation, for equal mass binaries, as a function of 
redshifted chirp mass. We can see that at 10 Gpcs, we should be able to see sources with redshifted chirp masses of 
10 s Mq, both with and without the extra corrections. The HHCs give median improvements of factors of 26 in chirp 
mass, 194 in reduced mass, 17 in time to coalescence and 12 in sky resolution. However, we should adapt a cautionary 
air about these sources. While they are detectable, it does not look like we will be say much about parameter 
estimation. Two of the most important quantities for astronomical observations are sky position and distance. For 
these particular sources, the distance to the source is unresolvable. The factor of 12 improvement in sky resolution 
minimizes the error from all-sky to 1 steradian. Again we should point out that these are median values, so there will 
be a small subset of the results which will give bigger improvements in the position and distance. But even then, it 
may not be good enough for astronomical purposes. 

So it is more realistic in this scenario to focus on the lower mass systems for parameter estimation. In this case 
we obtain more modest results, except in the case of the reduced mass. From the Monte Carlo we obtain HHC 
improvements in chirp mass of about 3.5 at 10 7 Mq down to approximately 1.7 at 10 4 M Q . For coalescence time we 
obtain improvements of 2.7 at high mass to 1 at low mass and an improvement in sky resolution of 5.3 at high mass 
to 1.3 at low mass. The smallest improvement comes for luminosity distance with 1.3 at 10 7 Mq to 1 at 10 4 M Q . 
However, the biggest improvement is for the reduced mass. Here we get a factor of 28 improvement at high mass, 
increasing to a factor of 73 at low mass. We will explain this fantastic result in the next section when we look at 
correlation breaking between the parameters. It is clear from the plot that the biggest gain from using the HHC 
waveforms is in the redshifted chirp mass range of 10 7 < M c /Mq < 5 x 10 5 . Below this value, the higher harmonics 
are now out of band at high frequencies. This leaves only the Hjl^ correction in band. From the cell where we 
display SNR, we can see that this first correction has the effect of subtracting SNR. The result that the HHC waveform 
has greater SNR at high mass, and lower SNR at low mass is well known. 

B. Unequal mass binaries. 

In Fig (|2|) we plot the same values for the mass ratio 10 scenario. In this case we are not able to detect binaries with 
a redshifted chirp mass of 10 8 -Mq, even with the inclusion of HHCs. For unequal mass binaries, the effect of the HHCs 
is somewhat less dramatic for the mass parameters. In fact we can see that we only have improvements in parameter 
estimation above 10 6 Af Q . Once again, between 10 7 M Q and 10 4 M Q , we have improvement factors of approximately 
3 to 1 for both chirp mass and reduced mass. However, for t c , A57 and D^, there is an obvious advantage to having 
the HHCs. While the improvement is not great (4 ~ 1), we do see a gain in the time to coalescence estimate. The 
improvements for luminosity distance are between 3.5 and 1.5, but the greatest improvement is for the sky resolution 
with median factors of 24 at 10 7 M Q and 1.5 at 1O 4 M . What is important in this case is that for both of these 
parameters, the improvement stretches to pretty much all lower masses. In terms of SNR, we see the same pattern as 
before. At higher masses the HHC waveforms have a slightly higher SNR, whereas the opposite is true at low mass. 

We should say a few words here about why the HHCs are having an effect. This was outlined in Refs 0, [f|, but 
we feel it is useful to reiterate it here. We have seen that our measured signal is a function of nine parameters. The 
three most important, and hence easiest to measure, parameters are {M c ,/z, t c }. We can see from Eqn Q that the 
waveform phase, and from uj = dQ/dt, that the frequency and all its derivatives are determined with these three 
parameters. To get good a good estimate for sky position we rely on two effects : first, the Doppler shift which is a 
function of {9, cf>} and the beam pattern functions which are a function of {9, 4>, ip}. For a system like the inspiral of 
a massive black hole binary, in effect each detector measures an amplitude and a phase. Because the waveforms are 
long-lived in the detector, we are able to measure the frequency and it's derivatives to high accuracy. This allows us 
to measure {M c , fj,, t c } quite accurately, but still leaves six unknowns {9, cj>, tp, i, (p c , \ii(Dl)}. 

In a two detector system with no HHCs, each detector measures a phase and an amplitude, giving in effect, four 
observables with six unknowns. However, the introduction of the HHCs rectifies this. In the equal mass case, we 
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FIG. 3: Main correlation breaking for equal mass binaries with a chirp mass of 10 7 Mq due to the inclusion of higher harmonic 
corrections. The restricted waveform correlations are given by the dark line and the HHC correlations are given by the lighter 
line. 



stated that all harmonics which are a function of dm = mi — mi are null. This means that if we include all harmonic 

(2) 

corrections up to H\ x , we still retain three harmonics, giving us three phases and three amplitudes in each detector. 
This gives us twelve observables in total for six unknowns. For the unequal mass case, we now end up with twenty-four 
observables for six unknowns. We can see the effect of having the extra observables by noting that the parameter 
estimation is always better in the unequal mass case. 



V. CORRELATION BREAKING DUE TO THE CORRECTED WAVEFORMS. 



In the previous section we demonstrated that the HHCs lead to smaller predicted errors in the estimation of 
parameters, especially in the higher mass range. While we would expect as improvement due to the fact that the 
HHC waveforms extend to higher frequencies, this still does not answer the question of what exactly is causing the 
improvement. To try and give a more definitive reason, we look at correlation breaking between parameters. 

We stated earlier that the inverse of the FIM, Fy, is the variance-covariance matrix S y . Rather than working with 
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FIG. 4: Main correlation breaking for unequal mass binaries with a chirp mass of 10 7 Mq due to the inclusion of higher harmonic 
corrections. The restricted waveform correlations are given by the dark line and the HHC correlations are given by the lighter 
line. 



this, it is more informative to work with the matrix of standard deviations and correlation coefficients defined by 

!V& i = j 

(26) 

In this matrix, the diagonal elements range between [—1, 1], where 1 is perfect correlation, -1 is perfect anticorrelation 
and represents no correlation. For supermassive black holes there is some degree of correlation between most of the 
parameters. However, in most cases the correlations are mild. As we are concerned with any correlation breaking 
between the strongly correlated parameters, we will focus on parameter correlations with an absolute value of > 0.5. 
As we have seen that the HHCs have a greater effect on the higher mass binaries, for this exercise we will only concern 
ourselves with the M c = 10 7 Mq case. Again, we will treat the equal and unequal mass cases separately. 



A. Equal Mass Binaries. 

In Fig we plot the strongly correlated parameter breaking for the case of equal mass binaries. The dark lines 
represent the restricted waveforms, while the light lines denote the HHC waveforms. As was shown in Ref [Til Il2l. |25| . 
there is a strong correlation between {M c , /i, t c } in the restricted PN case. What was not clear in these previous 
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works is the fact that these three parameters are also quite highly correlated with the azimuthal sky angle <p. On the 
other hand, it turns out that cos 9 is highly correlated with the three extrinsic parameters {cos l, Dl, </? c }. Finally in 
the restricted case, one of the other parameters of interest is the polarization angle ip which is correlated mostly with 
{cos L, (fi c } 

If we now focus on the light curves, we can see the effect of introducing the higher harmonic. The most surprising 
thing we can immediately see is that for {M c ,/j,}, rather than breaking the correlation, the HHCs actual make the 
reduced mass perfectly correlated with the chirp mass. In fact, the Fisher elements for these two parameters are 
numerically equivalent in the equal mass case. As the chirp mass is already well resolved, even in the restricted 
waveform case, this perfect correlation explains the huge increase in precision for the reduced mass in Fig {l}. Going 
back to the top two cells in this figure, we can see that the error curves for both parameters are almost identical. 

So while there is an increase in correlation for {M c ,/j,}, we can see that there is a huge correlation breaking for all 
the other intrinsic parameters upon introduction of the HHCs. We can see that both M c and /i decouple from t c , with 
median absolute correlations dropping from (0.46, 0.78) to (0.18, 0.18) respectively. In terms of the improvement in 
sky resolution, it looks like there is a combined effect from the higher harmonics. 

The first is that the azimuthal sky parameter (p decouples from {t c , /j,} with correlations dropping from (0.55, 
0.56) to (0.24, 0.08) respectively. Secondly, the polar angle cos9 decouples from {cos l, Lp c }. Here the correlations 
decrease from (0.66, 0.58) to (0.49, 0.48). Finally we see the same thing happening with the polarization angle ip with 
correlation decreasing from (0.55, 0.81) to (0.32, 0.66) respectively. We believe these various breaking of correlation 
between {cos (9, (p, ip} and the other parameters are the main reason for the increase in sky resolution. 

B. Unequal Mass Binaries. 

In Fig Q we plot the same thing for the unequal mass case. Here we find a similar story to the equal mass 
case. The main difference is that in this case we do observe some slight correlation breaking between M c and fj, 
with the correlation reducing from 0.94 to 0.87. Once again, the main intrinsic correlation is between {M c , /i, t c , </>}. 
Firstly, both M c and /z decouple from t c with correlations dropping from (0.83, 0.96) to (0.24, 0.22) respectively. This 
is accompanied by correlation breaking between {M c ,/i,i c } and <p with respective correlations of (0.83, 0.97, 0.98) 
dropping to (0.4, 0.68, 0.54). 

Again we see a drop in correlation between cos 9 and {cos t, Dl, f c } with values changing from (0.53, 0.8, 0.46) to 
(0.19, 0.42, 0.19). And finally, as in the equal mass case, we also observe some correlation breaking between ip and 
{cost, ip c } with correlations being reduced from (0.52, 0.79) to (0.28, 0.68) respectively. 

VI. A ONE VERSUS TWO CHANNEL LISA, WITH HIGHER HARMONIC CORRECTIONS. 

The final issue we look at is the advantage of using HHC waveforms if there is a problem with LISA and we lose two 
detector links. While the goal is always to have six links working, the question stands : can we still do astronomy with 
a broken detector? In the following text we will refer to one channel restricted results as Xq, two channel restricted 
as AEq, one channel HHC as X-i and two channel HHC results as AE%. There are two results we can present here, 
before making distinctions between equal and unequal masses. The first is that a full LISA is always going to perform 
better, and the second is that there is in general an improvement in SNR of ~ y/2 in going from a one to a two channel 
LISA regardless of the waveform type. 

A. Equal Mass Binaries. 

In Fig (O, we repeat the results of Fig (JXJ) , this time with the inclusion of parameter estimation just using one LISA 
channel. We can see that in all cases, losing one LISA channel and just having the Xq output would seriously hamper 
our ability to carry out parameter estimation for all parameters. For the M c = 10 8 Mq case, while the signal would 
be just about detectable over a galaxy with a median SNR of 12.6, we could say nothing about {M c , fi, Afi, Dl} and 
would have an enormous error in time to coalescence prediction. As we go to the M c — 10 7 M case, the gap between 
the Xq and AEq results have narrowed for {M c , fj,,t c } with a factor of ~ 2 in error. However for sky resolution and 
luminosity distance there is a loss of information of factors of 27 and 4.5 respectively. This is again due to the fact 
that in one channel, we would only measure one phase and one amplitude, giving two observables for the six unknowns 
{9,(P^,l, Vc MDl)}- 

With the introduction of the HHC waveforms, it is clear, as we can see from the plots, that for a chirp mass of 
> 2 x 10 6 M Q we can do as well, if not better, using the X2 output as we can using the AEq channels. Again at 
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FIG. 5: A comparison of parameter extraction for equal mass binaries using restricted one channel (dotted lines), restricted 
two channel (solid line) and higher harmonic corrected one channel (dash-dot line) and two channel (dashed line) waveforms, 
as a function of redshifted chirp mass and at a distance of 10 Gpc. The values quoted are median values from a 2 x 10 4 point 
Monte Carlo simulation. 



M c = 10 8 Mq, the sky position and luminosity distance are unresolvable using the X2 channel only. However, for 
chirp mass, reduced mass and coalescence time we actually improve the error estimate by factors of almost 10, 57 and 
4 respectively over the AEq channels. At around M c = 10 7 M Q we see the real power of including the HHCs. Except 
for luminosity distance, where there is a slight increase in error, the performance of the X2 channel equals or exceeds 
the AEq channels. In fact, we see improvement factors of 2 in chirp mass and 16 in reduced mass when we use just 
the X2 channel. We can justify this increase in performance by examining the observable count. We stated earlier 
that the AEq combination measures two phases and two amplitudes, giving four observables for six unknowns. Once 
we introduce the HHCs in the equal mass case, we now have three amplitudes and three phases, giving us the six 
observables we need for the six unknowns. 



B. Unequal Mass Binaries. 

We present the same results for the unequal mass case in Fig ©. As in the equal mass case, we see some similar 
patterns. In all cases an Xq channel alone will seriously effect our ability to do GW astronomy, with the main effect 
coming in sky resolution and luminosity distance estimations. We can see that the errors for these two parameters 
changes by at least an order of magnitude, while there is a loss of almost an order of magnitude for the mass 
parameters. For the X2 channel case, we again see that at M c > 2 x 10 6 M the one channel with HHCs does as 
well if not better than the AEq channels. As is expected, the improvement in the masses is not as pronounced in the 
unequal mass case, but we still see improvements at the highest mass. The major result in this case is that an X2 
channel will outperform the AEq cases down to a redshifted chirp mass of ~ 10 5 M & . 

It is interesting to examine the effects of the HHCs in the one detector case in a deeper manner. To do so, we will 
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FIG. 6: A comparison of parameter extraction for binaries with a mass ratio of 10 using restricted one channel (dotted lines), 
restricted two channel (solid line) and higher harmonic corrected one channel (dash-dot line) and two channel (dashed line) 
waveforms, as a function of redshifted chirp mass and at a distance of 10 Gpc. The values quoted are median values from a 
2 x 10 4 point Monte Carlo simulation. 



focus on the the case of a system with a redshifted chirp mass of 10 7 M© and a mass ratio of 10. As usual, the source 
is at 10 Gpc and is coalescing just inside the observation time. The first thing we can look at is the observables count. 
In this case, the extra harmonics allow us to measure six phases and six amplitudes giving us twelve observables for 
six unknowns. In fact, it is now clear that just adding the terms are enough to outperform the AEq channels, 

as with this model we would have three phases and three amplitudes, giving us the six observables needed for the six 
unknowns. We should note however, that for the equal mass system, we would need to go to higher harmonic orders 
due to the null harmonic terms. 

In Table [J we present the median errors for a Monte Carlo simulation of 2 x 10 4 systems. We can see from having 
the extra observables, how the X2 channel fares against the Xq channel, and more importantly, the AEq channels. 
As long as we include the HHCs, a one detector LISA would still measure {M c , fi, t c , cos 8, ip, \u(Dl)} factors of 
(1.5,1.7,1.9,1.5,1.3,1.4) times better than a two detector LISA without. We should point out here that one of the most 
important improvement comes in the measurement of the polarization angle ip. In Fig J7]) we plot histograms for the 
error in the estimation of ip with and without HHCs. We can see that having the extra HHCs allows us to make 
an instantaneous measurement of ip which is important firstly, in determining the binaries principle polarization axis 
around the line of sight, and secondly, in the overall amplitude of the detector response. 

The final thing we will investigate is the correlations between the parameters. It is interesting to trace the evolution 
in correlations to see the effects of the inclusion of HHCs and then the inclusion of a second channel. In Fig (jHJ), we 
display some of the main correlation breaking with a single detector. In Table [TTI we tabulate the correlations in both 
the one-detector case with and without HHCs, and in the two detector case with HHCs. As displayed in previous 
sections, the parameters {M c , /j,, t c , <p} are all highly correlated with each other. However, the introduction of the 
HHCs again cause a decoupling between these parameters. As well as that we also notice a decoupling between cos 8 
and Dl and importantly, between <p and ip. It would seem that the decoupling of <fi with four other parameters is 
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FIG. 7: Histograms for the error in the polarization angle ip for restricted (dark) and corrected (light) waveforms using a single 
LISA detector, for systems with a redshifted chirpmass of 10 7 Mq, a mass ratio of 10, at a distance of 10 Gpc. 
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FIG. 8: An examination of parameter correlations for the case of a single LISA detector for systems with a redshifted chirpmass 
of 10 7 M Q , a mass ratio of 10, at a distance of 10 Gpc. The dark lines denote the restricted waveforms, while the light lines 
signify the corrected waveforms. 



13 





ln(M c ) 


ln( M ) 




HQ 




cos 8 


4> 


cos I 


HDl) 


i> 






8.3 x 10" 4 


7.26 x 10" 


-3 


4.1 x 10" 


5 


0.269 


2.089 


0.118 


0.412 


0.429 


2.06 


AX 2 


3.2 x 1(T 4 


2.23 x 10" 


-3 


1.2 x 10" 


-5 


0.029 


0.315 


0.078 


0.095 


0.077 


0.81 


AAE 


4.73 x 10~ 4 


3.83 x 10" 


-3 


2.3 x 10" 


-5 


0.043 


0.06 


0.048 


0.12 


0.11 


0.3 


AAE 2 


1.65 x 1(T 4 


1.12 x 10" 


-3 


6.1 x 10" 


-6 


0.0203 


0.037 


0.009 


0.028 


0.029 


0.169 


X0/X2 


2.59 


3.26 




3.3 




9.44 


6.63 


1.51 


4.35 


5.57 


2.54 


AE0/AE2 


2.89 


3.42 




3.8 




2.1 


1.62 


5.33 


4.3 


3.8 


1.78 


Xq/AEo 


1.76 


1.89 




1.74 




6.26 


34.8 


2.46 


3.4 


3.9 


6.85 


X2/AE2 


1.8 


1.99 




2.0 




1.4 


8.5 


8.67 


3.34 


2.6 


4.79 



TABLE I: Median errors for a source with a redshifted chirpmass of 10 7 Mq and a mass ratio of 10 at a distance of 10 Gpc. 
For each parameter we show the one-channel LISA correlations without HHCs, Xo, one-channel LISA correlations with HHCs, 
X2, two-channel LISA correlations without HHCs, AEo, and two-channel LISA correlations with HHCs, AE2, errors. The last 
four rows show the ratios of the median error values for various channel combinations. The errors for the angular variables are 
in radians. 
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TABLE II: Absolute median values of the correlation matrix for a source with a redshifted chirpmass of 10 7 Mq and a mass ratio 
of 10 at a distance of 10 Gpc. For each parameter we show the one-channel LISA correlations without HHCs, Xo, one-channel 
LISA correlations with HHCs, X2, and two-channel LISA correlations with HHCs, AE2- 

responsible for the improvement in sky resolution at high masses. 

When we move to the two detector case with HHCs, we can see that the extra channel now serves to refine the 
parameter estimation. This is due to a slight further decrease in some of the parameter correlations. While we have 
not shown the results here, the introduction of a second channel only causes a slight decrease in the breaking of 
correlations. It is clear that it is the HHCs that are having the main effect. 
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VII. CONCLUSION 

In this work we have looked at the effect of including higher harmonic corrections to the restricted PN waveforms in 
the LISA context. By carrying out an extensive Monte Carlo simulation for various redshifted chirp masses and mass 
ratios, we have tried to cover as many scenarios as possible. As was already known, the corrected waveforms bring 
previously invisible sources into the LISA bandwidth. We have shown however, that in general while we will be able 
to detect these sources, we will not be able to say anything useful for any electromagnetic follow-up. For equal mass 
binaries, due to a combination of correlation breaking and the reduced mass becoming perfectly correlated with the 
chirp mass, there is a significant improvement in mass estimation using the corrected waveforms. We found that while 
there are systems where improvements in sky resolution and luminosity distance are huge, the median improvements 
are not as dramatic as some individual sources. The same is true in the unequal mass case, although here there is an 
order of magnitude improvement in sky resolution. We also showed that the overall effect of the HHCs is providing 
enough observables to account for the number of unknown parameters. 

We finally looked at the more interesting case where the LISA output was reduced to one channel. We showed 
here that for supermassive and some massive systems, we can actually do better with a one channel LISA with HHCs 
than we can with a full LISA with no harmonic corrections. As the HHCs again provide enough observables to solve 
for the unknown parameters, this allows us to improve mass measurement, sky resolution and make instantaneous 
measurements of the polarization ip. While a second channel always improves the parameter extraction, it is clear 
that the HHCs are the main source of correlation breaking and improvement in parameter estimation. This is an 
important result as it shows that while we would always like a full LISA, not all would be lost if we were forced to 
work with a single channel detector. 
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